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Abstract 

The threshold behavior of the master amplitudes for two loop sunrise self-mass 
graph is studied by solving the system of differential equations, which they satisfy. 
The expansion at the threshold of the master amplitudes is obtained analytically 
for arbitrary masses. 



PACS ll.lO.-z Field theory 

PACS ll.lO.Kk Field theories in dimensions other than four 
PACS ll.15.Bt General properties of perturbation theory 
PACS 12.20.Ds Specific calculations 
PACS 12.38.Bx Perturbative calculations 



1 Introduction. 



The sunrise graph (also known as sunset or London transport diagram) appears naturally, 
as a consequence of tensorial reduction, in several higher order calculations in gauge 
theories. Due to the presence of heavy quarks, vector bosons and Higgs particles all the 
internal lines may carry different masses, so that sunrise amplitudes depend in general on 
three different internal masses rrii, i = 1,2, 3, besides the external scalar variable p^, if 
is the external momentum (in n-dimensional Euclidean space). 

For a proper understanding of their analytical behaviour, as well as for a check of the 
numerical calculations, it is convenient to know the amplitudes off-shell, and also around 
some particular values of p^, such as p"^ = 0, p^ = oo, = —(mi + m2 — ms)^ (one of the 
pseudothresholds) and p"^ = —{mi + m2 + ms)^ (the threshold). 

This paper is devoted to the analytic evaluation of the coefficients of the expansion of 
the sunrise amplitudes in p^, at the threshold value — = (mi + m2 + ms)^ = sq . The 
approach relies on the exploitation of the information contained in the linear system of first 
order differential equations in p"^, which is known to be satisfied by the sunrise amplitudes 
themselves It is to be noted that all the above points { p"^ = 0, oo , threshold and 
pseudothresholds) correspond to the Fuchsian points of the differential equations, which 
therefore emerge as a natural tool for their discussion. 

The analytic properties of Feynman diagrams at threshold and pseudothresholds are 
well known, see for example . The sunrise diagram, with different masses, has been in- 
vestigated in 1^ , while in [|] the values of the amplitudes at threshold and pseudothreshold 
were obtained. With the method established in and |^ , further, the expansion around 
threshold was obtained in 0, even if a complete analytical result was given there only for 
the case of equal masses. With the configuration space technique the expansion around 
threshold was investigated also in |^. The expansion around one of the pseudothresh- 
olds (the others are straightforwardly given by a cyclic permutation of the masses) was 
obtained in analytical form in 0. 

The threshold is a Fuchsian point of the system of equations obtained in In this 
case it is known |T0|] that the master amplitudes can be expanded around that point as 
a combination of several terms, each equal to a leading power x°'' times a power series in 
X, with X = (p^ + (mi + 1712 + ^3)^), where the exponents are in general real numbers. 

When the expansions are inserted in the differential equations, the equations become a 
set of algebraic equations in the ai and in the coefficients of the expansions; the obtained 
algebraic equations can then be solved recursively, for arbitrary value of the dimension 
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n, once the initial conditions ( i.e. in the case considered here the values of the sunrise 
amplitudes at the threshold) are given. As those initial values are in turn functions of 
the masses, we find in our approach that the initial values themselves satisfy a system of 
linear differential equations in the masses. We expand in the dimension n around n = 4 
the equations in the masses and solve them explicitly up to the finite part in (ra — 4). This 
result presented in the next section for n = 4 and = —{mi +1712 + 7713)^ is in agreement 
with the literature @]. 

Once the initial conditions are obtained, we look at the recursive solution of the 
algebraic equations for the coefficients of the expansion in (p^ + sq) (our results are given 
in the Section 3 ). It turns out that the formula, expressing the second order coefficient of 
the [p^ + So) in the first master integral expansion in terms of the zeroth order values of 
the other master integrals, involves the coefficient l/{n — 4). Therefore the finite part at 
{n — 4) of the second order term in the {p^ + Sq) expansion involve the first order terms in 
(n — 4) of the zeroth order terms in (p^ + Sq). We evaluate those first order terms in {n — 4) 
solving a differential equation, which they satisfy, obtained also in Section 3. Fortunately, 
the formulae expressing the higher order coefficients of the (p^ + sq) expansion involve 
coefficients like — 5), l/(n — 6) etc., which are finite at n = 4, and can be used without 
further problems for evaluating those higher order terms. Our results agree in the equal 
mass case with those obtained in (we take a numerical constant from the comparison), 
while for not equal masses the result is given in an analytical form for the first time. 

2 The threshold values of the master amplitudes 

It is known that the two-loop sunrise self- mass graph with arbitrary masses mi, 7712, 
has four independent master amplitudes [jl2[, which will be referred to, as in [|l|, by 



where n is the continuous number of dimensions, mi,i = 1,2,3 the three masses and 
the external ra-momentum. FQ{n,m\,m\,m\,p'^) is the scalar amplitude 
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The values of the master integrals at the threshold can be obtained in a way similar to 
the way followed for the values at the pseudothreshold ^ and we will discuss here mostly 
those points of the derivation which are different from the pseudothreshold case. If the 
values at the threshold are written as 



ga{n,mi,m2,m3) = Fa(n,ml,ml,ml,p^ = -{mi+rn2 + m3)'^'j ,a = 0,l,2,3 (4) 
we find, expanding around n = 4, 



Gain, mi, 1712,1713) 



+g^^\mi,m2,m3)+0{n-A) a = 0,1, 2, 3 



(n - 4)' 



jGi ^\mi, 1712, 1713) + 
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^\mi,i7i2,m3) 
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The coefficient C (n) is a function of n only and at n = 4 has the following expansion 
C{n) = (2v/^)' 

= 1 - (n - 4) 
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+ (n-4)' 



\og\2^) + + 7s log(2v/^) 



+ 0((n-4f) , (6) 

where 7^; is the Euler-Mascheroni constant. 

From [0], where the singular parts in {n — 4) of Fa{n,mi,ml,ml,p'^) are given for 
arbitrary values of p"^, we have 



^\mi, 1712, 1713 

,(-1) 



{ml + m\ + ml) 



g^ '''{mi,m2,m3) = -^{mi + m2 + ms)^ + -^{mj + ml + ml) 

62 lb 



m\ log(m^) + m\ log(m2) + m\ log(m3) 

1 1 
"16 ^ 8 



g\ '^\mi,m2,m3) = - , g\ ^\mi,m2,m3) = -— + -\og{ml) , i = l,2,3. (7) 



To obtain the other coefficients of Eq.(|^), as in |^, our starting point is the sys- 
tem of four linear differential equations in p^, satisfied by the four master amplitudes 
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Fa{n, ml, 1712, 171^, p^) , a = 0, 1,2,3, which are not rewritten here for short. By imposing 
the condition Eq. (Q) to p"^ the system becomes equivalent to a system of three hnear dif- 
ferential equations in m^, which in turn can be written as a single third order differential 
equation in for the function QQ^\mi, 1122,7123). We do not repeat here the details of 
the derivation, which is very similar to the derivation leading to the analogous Eq.(17) 
of 1^ (with the important difference of the substitution —>■ —m^ ). Its solution can 
also be obtained again in an analogous way, for positive values of mi,m2 and m^; it con- 
tains three unknown constants of integration (indeed functions of the masses mi and 7712) 
Ciijni, 7712), i = 1, 2, 3, to be fixed by the initial conditions. 

The initial conditions have to be imposed however in a way different with respect to 
the pseudothreshold. The point = 0, which is known and was used to impose initial 
conditions at the pseudothreshold , corresponds in the case of the threshold to the value 
ms = — (mi+m2), which is outside the validity of the present solution for ^g^"* (mi, m2, m.3), 
where positive values of the masses are assumed. However we know already the value of the 
function at the pseudothreshold ^ and for ms = the threshold and pseudothreshold 
are identical. That allows us to choose that point for the initial conditions. Imposing 
that, at 7713 = 0, G'o°'*(mi, m2, ms), given by the Eq.(25) of 0, and its first derivative 
with respect to ma are equal to the function QQ'\mi,m2, m^) and to its first ms derivative 
respectively, one gets two relations between the three functions Cj(mi, m2), i = 1,2, 3. The 
second derivative at that point is infinite and does not give an additional relation. After 
eliminating the two coefficients Cj(mi, m2), z = 2, 3 one gets^r(mi,m2,m3) as a function 
of Ci(mi, m2) only. The remaining constant can be fixed by the relations obtained by the 
permutation of the masses as we know that the Qo'^imi, m2, m.3) is symmetric under that 
exchange. Once all the constants are fixed the solution reads: 



gf\mi,m2,m2) - 
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+ g("^3 -m^- TTL^) {Ct{mi, 1712, 1713) + Ct{m2,mi,m3)) 
1 

~32 



log (m^) + 1712 log ("^2) + "^3 log (^3) + ('^i + "^2 ~ ""^3) log("^i) log("^2 



+ {ml — m\ + m\) log(m^) log(m3) + {—m\ + m\ + m\) log(m2) log(m3 
— mi(7mi — 27713) log(^i) ~ I7i2{7i7i2 — ^m^) log(m2) 
+ (2m;^ + 2mim2 + 2m2 — Smg) log(m3 
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mim2m3{mi + m2 + ma) 



^^(m^ + m2 + ml) + ^(mim2 + mims + 71227123) 



where 
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777i7773 
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The above result is in agreement with To check it one has to expand in 77 around 
77 = 4 all the factors which depend on 77, such as our C(77) ( Eq.(^) and the 77-dependent 
functions of 0. 

The functions ^1(77, 7/71, 7772, 7773) i = 1, 2, 3 can be easily found in an analogous way 
as for the pseudothreshold and the solutions read 

^1(77,7771,7772,7773) = 

lid d 
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Q2{n,mi,m2,m3) = 

lid d 
2 ^ + "^2) -^zr^oin,mi,m2,m3) + (7772 + 37773) ^—^o('^, "71, 7772, 7773) 



57775 
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where T(n, m^) is defined by 
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• (14) 



The other two functions, ^1°^ (7771, 7772, 7773) and ^2°'' (^15 m2, 7773), can be easily obtained 
by a permutation of the masses and we do not report them here. Again the results are in 
agreement with |Q. 



6 



3 The expansion at the threshold. 



The expansions at the threshold of the master amphtudes can be found with the help 
of the system of equations obtained in []T| (Eqs.(5,7)). As the threshold is a singular 
point the expansions consist of the regular parts plus the singular ones with fractional 
exponents. By inserting the expansion into the system of equations one finds that the 
fractional powers are fixed and the expansion reads 

CX) 

where n is the continuous dimension, a = 0, 1, 2, 3; j(0) = 2, j(l) = j(2) = j(3) = 3 and 

X = p'^ + (mi + 1712 + f^sY (16) 
is the expansion parameter. Of course in this new notation 

7Y'^°'°)(n,mi,m2,m3) = Qain, 1711,1712,1713) (17) 

is already presented in Eq.(^. 

From Eq.(^ it follows that the 7Y''"'*^(n, mi, m2, 7713) and 'H[°''^\n, rrii, m2, m^) are re- 
lated, so that having H^^''^\n, mi, m2, ^3) and H^^'^\n, mi, m2, m.3) one can calculate the 
others for a = 1, 2, 3. 

From the discussed system of equations follows also that the following differential 
equations hold for Hf''^^ {n, mi, m2, m3) 



n, 7711,7712 , 7713) 



9m 1 

dHf'^^Ti, mi, m2, m3) 
9m2 

dT-Cf'^^n, mi, m2, m3) 
c?m3 



71 — 3 



+ 



5 — 3n 



2mi 2 (mi + m2 + m3) 



71 — 3 



+ 



5 — 3n 



2m2 2 (mi + m2 + m3) 



71 — 3 



+ 



5 — 3n 



2m3 2 (mi + m2 + ms) 



Hf'°\7l, 7711, 7712, 7713) 

7Yf'°^(n,mi,m2,m3) 
7Yf'°)(n, mi,m2,m3) 



That means that 1-Lf'^\7i,7ni,7n2,7n3) can be written as 
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[mim2m3 



n — 3 
2 



(mi + 1712 + "msY 



(19) 



where C-h {n) is a function of n only, to be discussed later. As usual, the coefficients 
7i^°'*)(n, mi, m2, ma), i = 1, ■ ■ ■ can be found by solving a system of (in this case four) 
linear equations. We report here for brevity only one of the coefficients 



7^^°'^)(n, mi,m2,m3 
1 

16(mi + m2 + ma) 



3 (3n - 5) 

mi + m2 + m3 



, , 1 1 1 

in-3 — + — + — 
\mi m2 m3 
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The singular part of the expansion in Eq. ([l5D can be compared with the results of 0, 
as it corresponds, for arbitrary masses and in the language of that reference, to the (p-p) 
regions contribution completely expanded at the threshold. Our results Eq.(|l^), Eq.(|19D 
and Eq.(^) agree analytically with those presented in Eq.(5), Eq.(49), Eq.(50) of 0, 
provided that the value given in Eq.(49) of is taken with a minus sign, in agreement 
with its equal mass limit, Eq.(51) of 0. 

In the regular part of the expansion all the coefficients can be found provided we know 
?-^("'°)(n, mi, m2, ms) = Qai^, mi, m2, m^), for a = 0, 1, 2, 3, already given in Eq.(^. This 
is only partly true as we know 7Y''"''^''(n, mi, m2, ma) expanded around n = 4 up to the 
constant term only. The next to lowest term in the x-expansion reads 



H^'^'^^ {n, mi, m2, m^) 



1 



[n - 3)go{n, mi, 1712,1713) 



(mi + m2 + ma)^ 

+mlgi{n, mi, m2, ms) +mlQ2{n, mi, m2,m3) +mlg3{n, mi, 7712,7173) . (21) 

The Eq.(|21|) allows us to find the expansion of Ti^^'^^n, mi,m2, ma) around = 4 up 
to the constant term as we know the expansion of Qa{n,mi,m2,m3) up to the constant 
terms. It can be written as 



l-6^'^\n,mi,m2,m3) = C^(n) 



(0) 
(0,1) 



0(n-4) 
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«(o,i) - 32 ' 



(23) 



and 
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- 4(mi + m2) (mlCt{mi,in2, iris) + m\Ct{m2, mi, ms) 
+ (^m^ + mim2 + (^mi log(m^) + m2 log 

— {m\ + 2m\m2 + 2mim2 + log 



[rrin 
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32(mi + m2 + ma)^ 



4m^£t(mi, m2, ms) + 4m2A(m2, mi, ma) — 2(m^ + mim2 + m\) 



— mi(3mi + 2m2) log(m^) — m2(2mi + 3m2) log(m2) + (3m^ + 4mim2 + 3m2) log(m3 
mi + m2 



16(mi + m2 + ma) 32 



+ ^ log(m3) - 



128 



(24) 



In the equal mass case (mi = m2 = m^) the above result is in agreement, up to terms 
0(n — 4), with , while the analytical result for all masses different was not previously 
known in the literature. 

In order to obtain the first order term in the x-expansion for the other master ampli- 
tudes, we choose to proceed with the expansion of Fo(n, m\, m\,p^). It is necessary to 
know terms of the order of -Fo(^; ""^i? "^2; "^ii P^) find terms of the order x of the other 
master amplitudes using the relations of Eq.(^). For the expression of H^^''^^ {n, mi, m2, m^) 
we find 



7i^°'^^(n, mi,m2,m3) 



32(n — 4) (mi + m2 + ma)^ 
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(2mi + m2 + ni'i)Qi{n, mi, m2, ma) + (mi + 2m2 + m3)Q2{n, mi, 1712, ms) 



+ (mi + m2 + 2m3)Qs{n, mi, m2, m^) 



(25) 



We do not write here explicitly other terms, which are known, containing up to triple 
pole in the (n — 4) expansion. In fact as one can see later (Eq.(^)) the triple and double 
pole terms cancel in the expansion. The presence of the explicitly shown term ~ is a 
kind of an obstacle analogous to the one encountered at the pseudothreshold expansion 
0. In principle it requires the knowledge of the expansion of the combination of the 
master integrals written in square brackets of Eq.(^) up to the term ~ (n — 4) included. 
We could solve that problem in an analogous way as in the pseudothreshold expansion , 
but to show how far we can get just relying on the differential equations we will proceed 
differently. Defining the above combination 



gc{n,mi,m2,m3) 



(2mi + m2 + m3)^i(?T,, mi, m2, ms) 



+ (mi + 2m2 + m3)Q2{n, mi,m2, m^) + (mi + m2 + 2m3)^3(n, mi, m2, m^) 
we calculate the m^ derivative of it 



(26) 



dQc{n,mi,m2,m3) 
dm^ 



{n - 3)^c(n, mi, m2, m^) (-^ + — — 

^2m3 2(mi + rr 

+gc{n,mi,m2,m3)( — — 

\ mo 



1^2 + m^) mi + m2 
1 2 

m3 mi + m2 + m3 mi + m2 + 2m3 

/ 2mi+m2 3mi+m2 

+ {n - 4)gi(n,mi,m2,m3) -1 \ ■ — 

^ m3 mi + m2 + 2m3 

mi + 2m2 mi + 3m2 



+ 2m2 



+ {n - 4)^2(n, mi, m2, m3) (^-1 - 
+ ^'!~^^' T{n,ml)T{n,ml){l 



ms 



mi + m2 + 2m3 
mi 



4m2m3 

^^T(n,m?)r(n,m^)(l 



m2 

2m3 2(mi+m2 + m3) 
mi m2 



2m3 2(mi+m2 + m3^ 



+ 2^2 ^(^' mDT{n, ml) ( 



2,/ mi+m2 



4mim2 



+ 



mi + m2 



2m3 2(mi + m2 + ma) 



(27) 



10 



When expanding around n = 4, we find (due to the factor (n — 4) in front of 
Qi{n, mi, 1712, TTT's) and Q2{'n, mi, m2, m^) ) that the term of Qd"!^, mi,m2, m^) proportional 
to (n — 4), denoted as G^^\mi,m2,ms), 



Gc{n,mi,m2,mz) =C'^{n) h (n - 4)G^^)(mi, m2, ms) H , (28) 

fulfills the following first order differential equation 

9G'W(mi,m2,m3) 1 n). , / 1 1 n 
^ = 7;Gy{mi,m2,m3)[ j 

(7777.3 2 ^7773 777 1 + 7772 + 7773^ 

,1^(0). /27771 + 7772 777i N 

2 ^ 7713 7711+7712 + ^3^ 

,1^(0). /27712 + 7711 7712 \ 

+ o^2 (^l,^2,m3)( + ■ ■ ) 

2 ^ 7773 777i + 7772 + ^3 ^ 

1 (0). / 3(7771 + 7772) 777i + 7772 \ 
+ 0^3 ' mi, 7772, 7773 (4 + ^ ^ ) 

2 ^ 7773 777i + 7772 + TTls ' 

^1 3/ o I 27771 + 7772 777i x 
— log (777i) —2 H I 



384 ^ 7773 777i + 7772 + ^,3 

777i + 27772 "72 



i-V(777^)(-2 + 



384 ^ 777.3 m.1 + 777.2 + m.3 

1 1 3/ 2n/ . 7771 + 7772 777 1 + 7772 \ 
— ^(7773) -4 H ^ ) 



384 ^ 7773 7771 + 7772 + 7773' 

1 1 / 2M f 2sf, / 2^ , , / 2n\/"^1+"^2 777 1 + 7772 N 

— log(777 J log(7772) ( log(7r7 J + log(7r72) ) I ■ ■ ) 

/.X ^ / ^ m-y m-i -+- mn -+- tdo ^ 



128 ' ' ' ' ^'^^ 7773 7771 + 7772 + 7773' 

^ l0g(777?) l0g(777^) (log(777?) + log(777^)) (-2 + ^ + ) 



128 ^ 7ri3 7771 + 7772 + ^3' 

-T^ log(7r7^) log(7r7^) (log(7r7^) + log(777^)) (-2 + ^ + ) , (29) 

128 ^ 7773 7771 + 7772 + 7773^ 

instead of being a part of a system of three differential equations, as in the case of 
arbitrary ti. 

We can find a solution of that equation in a relatively simple way, which consists 
mainly in integrating by parts, so we do not report here details of the derivation. The 
solution reads 



^(1)/ X (7772 - 7773) (27771 + 7772 +7773) x 

G), ' [mi, 7772, ^3) = ■ ■ r Ct[mi, 7772, "^3) 

4(777i + 7772 + ^.3) 
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1 



{nil -rris) {nil + 21712 + ms) ^ . . 

H 77 \ \ ^ Ct{m2, mi, ma) 

4(mi + m2 + ni3) 

log'^(m^) (2mi + m2 + m^) + log^(m2) (mi + 2m2 + ms) 



192 

+ log^(m3) (mi + m2 + 2m3) 



+^ log(mi) log(m2) (log(mi) + log(m2)) (mi + m2) 



64 
1 

64 



log(m^) log(m3) (log(m^) + log(m3) ) (mi + ms) 



+ ^ log (ma) log(m3) (log(m2) + log(m3)) (m2 + mg) 



32 
1 

+ 77 



ms 



(^log(m^) + log(m2)) + m2 (log(mi) + log(m3)) + mi (log{nil) + log(m3) 



log(m^) (— 3mi + m2 + ms) + log(m2) (mi — 3m2 + ms) 



+ log(m3) (mi + m2 - 3m3) + — (mi + m2 + ms) 
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16 



/ mim2m3 
mi + m2 + ms 



IT 



3 1 

- (log(m^) + log(m^) + log(m^)) + - log(mi + m2 + ms) 



- log(mi + m2) - log(mi + ms) - log(m2 + ms) 



Js(mi,m2,ms) - K 



(30) 



where the unknown function of mi and m2, which remains after integration, is reduced 
to the single constant K using the symmetry of G^^''(mi, m2, ms) under the interchange 
of all the masses. The actual value of K will be fixed later. Ct{nii,m2,ni3) is defined in 
Eq.(^ and the only non trivial integral left, Ts(mi, m2, ms), is 



Xs(mi, m2, ms) = Xs(mi, m2, ms) + Js(mi, mi, m2) - Xs(m2, mi, mi) 
Ts(mi,m2,ms) = v^mlm^ / rfms 



m^{mi + m2 + ms) 



log 


( m3 \ 




\mij 



+ 



mis + mi ms + m2 



t\log (^^^^^) [log(t - ti) - log(t - t2)] + log (^^^^^) [log(t + t2) - \0g{t + ti)] 

+ log(t - ti) [2 log(l - ti) - log(ti)] - log(t - t2) [2 log(l - t2) - log(t2)] 
- log(t + ti) [2 log(l + ti) - log(-ti)] + log(t + t2) [2 log(l + t2) - log(-t2)] 

t-ti\ ^. ( t-t2 



't-ti\ {t-t2 

-2 Li2 + 2 Lis ^ 



1-ti 



t2 



IA2 



tl 



Li2 



t2 
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+2 Li2 ( 



- 2 Lio 



-t + t2 



where t and ti 2 are defined as 



(31) 



t 



1,2 



+ 7712 + ms - 

a/777i + 77I2 + 7773 + ^/m^ ' 

7772 ~ f^i ± 2iy/mim2 

777i + 7772 



(32) 



Having G^-^^ (7771, 7772, ^773) we can find the 77 = 4 expansion of (77, 7771, 7772, 7773) up 
to the constant term. It reads 



(77, 7771, 7772, 7773) =C' 



77 



(77-4) 



(33) 



where 



TT ^^777177727773(7771 + 7772 + ^3) 

" 32 (7771 + m2 + 7773)4 



(34) 



and 



7i 



(0) 
(0,2) 



1 



8(7771 + 7772 + rUs)^ 
1 



+ 



777^(7771 + 2m2)Ct{mi, 7772, '"^'a) + m\{2mi + 7772)^(^^2, ■"^■1, ''^3) 



4(7771 + "^2) (mlCti'mi, 1712, mz) + m\Ct{m2, mi, 7773 



16(777i + 7r72 + 7^73)^ 

— 777i log(777^) ^777^ + 777 i7772 + 7772^ — 7772 log(7772) ^777^ + 777 i7772 + 7772^ 
+ log(7773) {rn\ + 2777^7772 + 2mim\ + 7772^ 

'' Am\Ct{mi,m2, m^) — Am\Ct{m2, mi,m^) 



32(777i + 7772 + 



+ 777i log(7ri^) (37r7i + 27772) + 1^2 log(7772) (2777i + 37772) 

— log(7773) ^3777^ + 4777i7772 + 37772^ + 2(777^ + 777 i7772 + 7772) 
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+ 



32(mi +1712 + m-i) 
1 



—mi \og{m\) — 1712 log(m2) + {mi + 7712) log(m3) — 2(mi + m2] 



64(mi + 1712 + 



T{ ^ mim2m-i{mi + m2 + rris) 
4(mi + 7712 + ^3)^ 



(^m1 + 17111712 + ITI2 



71 mim2in:i{mi + m2 + ins) 
4(mi + 1712 + ^3)^ 
171177121713(1711 + 1712 + ins) I 



(mi + 



64(mi + m2 + m^) 



lITlo 



- (log(mi) + log(m2) + log I 
-2(log(mi + m2) + log(mi + ms) + log(m2 + ms)) + log(mi + m2 + 1713) - 3 
+2X3(mi, m2, mg) -2K\ . 



(35) 



As we know only the expansion of TC^^''^\n,mi,m2,m3) around n = 4 and not its 
exact form for arbitrary n, to proceed, we write also the expansion in (n — 4) of the 
complete first master amplitude Fo(n, m^, m2, mg, p^). The singular and regular part are 
well distinct for continuous arbitrary n; in the expansion around n = 4, the singular part 
gives terms proportional to log(a;) plus other terms , without log(x), which mix with the 
terms coming from the regular part. Expanding Eq.(|T3|) in (n — 4) we find 



Fo{n,m^,m2,m^,p ) = C {n)< 



1 



[n — Ay 



1711,1712,1713) + ^—tQo ^\Tni, 1712, 1713) + Qo\mi, 1712, 1713) 

[n — 4) 



+ X 



+x" 



{11 -A) 



(0,1) 



n 



(0) 



(0,2) -^(0,2) (log(x) + |log 



171117121713 



+ b 



mim2m3[mi + 1712 + 7713 



(mi + m2 + ms)^/ / 

+ 0(n-4,x^) 



(mi + 7712 + m3)^ 

where the constant b comes from the expansion of C-j-i (n) around n = 4 



(36) 



Cn (n) = C\7i) 



TT 



(n-4) 32 



0(11-4) 



(37) 



14 



The requirement of the disappearance of the pole term in the coefficient of in Eg . (1361) 
and the knowledge of the pole term of 7i^°'^''(n, mi, m2, ms) fix the pole term in C-^ (n). 
That requirement fix also the absence of higher pole terms in C-j-i (n) . In the expansion 
Eq.(^) there are two still unknown constants: K in the expression of 7i(o\) and b from 
Cn {n). They cannot be fixed separately due to the n = 4 expansion, but only in the 
occurring combination (6 — K/32). As we are interested only in the n = 4 expansion 
of the master integrals, the knowledge of the combination of the constants appearing in 
Eq.(^) is enough to fix all the higher order terms in the threshold expansion. To fix the 
combination of the constants it is sufficient to know the term ~ of -fo('^) ""^i, ''^35 P^) 
for fixed values of the masses. Due to the factor in front of b and K, y^m^m^rn^, all 
masses have to be different from zero and so the simplest choice is the equal mass case 
(mi = m2 = ms = m). As the analytical result for equal mass case is available 0, we 
limited ourselves to its numerical check, which can be performed by using the dispersion 
relation representation of the master integral -Fo(^ = 4, m^, m^, m^,p^) 0. Its second 
derivative reads 



d'^Fo{n = 4, m^, m^, m^,p^) 



du— — Eo{u,m) , (3J 



where 



a(p2)2 J 8{U + p2)3 

9m^ 



EQ[u,m) = — / db . (39) 



u 

Am? 



At = — 9m^ the integral Eq.(^) is logarithmically divergent (in agreement with 
Eq.(il)), but the divergent part can be easily extracted by integrating by parts, providing 
us with the two needed leading terms in the threshold expansion of the integral 



9^Fo(n = 4, m^, m^, m^,p^) 7i"V3 j ^ Xg 



16 /*'°n^^Hs/^ + 



9m2 



fi4 2 — 

o^m"' \m^ 

1 , ^ 1 ^ !^ ^ i2§)^ _ ^ci, a) + OM . (41) 



864m2 m2 \ 972 432 648 / 1296m2 
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where Xe = + 9m^ corresponds to the equal mass limit mi = m2 = = m of x in 
Eq.(|l^) and the explicit value of the integral is taken from [0. The analytical result of 
0, in Eq.(^), agrees with the numerical value of the integral in Eq. (pOD , so we did not 
reevaluate analytically the integral. 

For extracting the combination of the constants (6 — K/32), we consider Eq.(^) in 
the equal mass limit mi = m2 = = m, and we compare its second derivative with 
Eq.(^), obtaining 



b = TT 

32 



■^ + |log(2))+ici.(| 



(42) 



to be used in Eq.(|36|) and in all the other n = 4 expansions of the master integrals. For 
completeness we report Eq.(|36D in the equal mass limit mi = m2 = ms = m in our 
notations 



Fo(n, m ,m ,m ,p ) = C (n) 



3m 



+ 



3m 



(3-41og( 



7r 



6 

V3 23 1 



n-4)2 32(ra-4) 

+ -TT^ (l5 + 121o"'-2 



m 



m^7r^/3 3m^ 



128 



m 



8 log^ (m^ 



32(n-4) 108 384 32 

2 



log(m^ 




- log (^) + - log(2) + - log(3) + -)- -Ch (-) - ^ 
4 ^\m^J 4 ^ 2 ^ 24/ 4 V3/ 8 



(43) 



The higher order terms in the expansion at threshold can be easily obtained alge- 
braically, but they are of interest only in case one wants to calculate the master ampli- 
tudes using the threshold expansion , which is not our purpose here. The expansion of 
all the master integrals up to terms ~ x is however necessary for the numerical solution 
of the system of equations obtained in [jl|. 

The threshold expansion of the remaining master integrals can be obtained using 
Eq.(D. One gets 



17/ 2 2 2 2 



[n 



4) 



mi, 1712, ms) 



(n-4) 



Till, 1712, m^j 
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(44) 



where Qj^\mi,m2,m^) are defined in Eq.(|^) and Eq.([T^) and 



'^(3,1) - 



IT 



^ mim2m^{mi + m2 + m^) 



(mi + m2) 



mim2m^{mi + m2 + m^) 
8m3(mi + m2 + m^y 

1 /mi + mgN 1 /mi + ma 

+ o log + o log 

2 \ 7712 J 2 V ms 



TT ( — log ( 

\2 Wmi + m2 + ms)^ 



1 /m2 + mg 
+ 1 + - log 

2 V mi 



8(mi + m2 + ms)"^ 
1 

16(mi + m2 + 771^)^ 
1 

16(mi + m2 + rn^y 



77ilCt{mi, 7712, rriz) + mlCt{m2, mi, ms) 



mi log(m^) + m2 log(m2) — (mi + m2) log(m3 



(45) 



The expansion of -Fi(n, m^, m^, m|,p^) and F2{7i, m^, m2, m|,p^) can be obtained by a 
permutation of the masses. 



4 Summary. 

In this paper we have presented the expansion of the 2-loop sunrise selfmass master 
amplitudes at the threshold = —(mi + m2 + 7713^. We define the expansion in Eq.(|15]); 
the values of the amplitudes at the threshold are given in Eq.(||), Eq.(0), Eq.(||) and 
Eq.(|l^). The first order terms in the threshold expansion of the master amplitudes at 
n = A are presented in Eq.(|2^) and Eq.(^), while only for the first master amplitude 
Fo[n = 4, m^, m^, m|,p^ = —(mi + m2 + ms)^) the second order term in the threshold 
expansion is presented in Eq. (|5BD . The higher order terms, which are not given explicitly 
here, can be easily found by solving recursively, at each order, a system of four algebraic 
linear equations. As said already in the expansion at the pseudothreshold cannot be 
simply deduced from the known expansion at the threshold and vice versa, even if at first 
sight they seem to be connected by the change of sign ms —fri^. In fact the analytic 
properties of the amplitudes are different at the two points: at the pseudothreshold the 
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sunrise amplitudes are regular, so that the solution of the system of equations [jll can be 
expanded as a single power series, while at the threshold the sunrise amplitudes develop 
a branch point and its expansion is indeed the sum of two series Eq.(|l^), |T^. 
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